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CM I Abstract: We introduce a class of finite tight frames called prime tight frames 

and prove some of their elementary properties. We show that finite tight frames 
can be written as a union of prime tight frames. In addition, we show that when 
they exist, equiangular tight frames are prime tight frames. We then examine two 
families of well-known tight frames vis-a-vis the notion of prime tight frames. 
More specifically, we characterize all prime harmonic tight frames and use this 
characterization to suggest effective analysis and synthesis computation strategies 
for such frames. Finally, we describe all prime frames constructed from the spec- 
tral tetris method. As a byproduct of this last result, we obtain a characterization 
of when the tetris construction works for redundancies below two. 



1 Introduction 

CM 

A frame for a finite dimensional vector space is a spanning set that is not necessarily a basis. 
Thus, a frame naturally gives rise to the decomposition and reconstruction of any vector in the 
rv~) ■ space. This simple fact has important consequences in many areas of science and engineering 
^ . in which frames now play an increasingly important role. We refer to (9l[IIl[T6l|20l|2Tl for an 
overview of frames and some of their applications. 
[ In this paper we are interested in the classification of the tight frames that remain tight after 

the deletion of frame vectors. We remind the reader that a tight frame for K N is a collection of 
vectors {(pi}f =l C K N satisfying 

^ 1 

_53; M\4 2 = £K*'W)I ' forallxGK^ 

i=l 

for some constant A > 0. Here K is either R or C. For (non-tight) frames, this problem is 
known as the erasure problem for frames that was first investigated by Goyal, Kovacevic, and 
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Kelner lfl"5ll . and by Casazza and Kovacevic [[8]|. In particular, the focus in these works was 
whether any set of a given size can be remove from a frame to still leave a frame. In this case, 
estimates for frame bounds after erasures were given. Recently, a more general framework for 
the erasure problem was proposed by Bodmann and Paulsen in |@]. 

As mentioned above, we classify tight frames that remain tight after some erasure. In the 
process, we define a new class of tight frames called prime tight frames as tight frames for 
which no proper sub-collection is a tight frame. We show in Section [3] that prime tight frames 
are fundamental building blocks for all tight frames in the sense that any tight frame can be 
written as a union of prime tight frames (Theorem 13.41) . This is in a way analog to the prime 
factorization of natural numbers! We will elaborate more on this later in this paper, and simply 
mention here as an example, that all harmonic tight frames, i.e., tight frames constructed from 
anMxM discrete Fourier matrix by keeping the first N rows, are prime whenever M is a prime 
integer. In Section 14.21 we will see that the reason for this is that every Mth roots of unity are 
primitive and therefore no proper subset of these roots add to zero. We refer to 11241 1271 for 
results on vanishing sums of roots of unity. 

To further illustrate and motivate our results, consider a tight frame obtained through the 
spectral tetris method which was recently introduced by Casazza, Fickus, Mixon, Wang, and 
Zhou [|6l|. This method leads to the construction of tight frames of M unit norm vectors in M. N 
for any given N,M £ N satisfying M > IN. One of the key features of this class of frames is 
that they are sparse in the sense that the coefficient vector of each frame element with respect to 
an orthonormal basis contains only few nonzero entries O. For N = 4 and M = 1 1, the spectral 
tetris construction yields the frame <I> = {<P;} 6 ; with / = {1,2,. ..,11}: 




where the frame vectors of <I> appears as columns in the matrix. Observe that the first two vec- 
tors are identical, which might be an undesirable for encoding because they lead to transform 
coefficients (/, (pi) carrying no new information lfl"5Tl . Write $ = $i U 3»2 where 3»i = {<Pi} ie j , 
h = {1,5,8,11}, and 4> 2 = {<Pi} i€l2 ,h = {2,3,4,6,7,9, 10}. Notice that 4> 2 is a tight frame 
with redundancy reduced to 7/4, and which can not be reduced further, while 4>i is an or- 
thonormal basis. Consequently, 

11 4 

/ = IT E </> <Pi) <Pi = IT E (/> e >) e ' + TT E if' ^ = Tlf + TT E if' 
j=l i'=l ieh ieh 

hence / = f £ ; - €/2 (/, (pi) (pi for every / £ R 4 . 

The reduced tight frame 4>2 does not have the artifact of $ with respect to repetition of 
vectors. Moreover, it is a tight frame with redundancy M/N less than two, something that is not 
possible using spectral tetris algorithms from Hence, this reduced tight frame possesses 

additional desirable properties as compared to the original tight frame, without loosing the spar- 
sity of the original frame <£. We shall see in Section |4~41 that all spectral tetris tight frames can 
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be decomposed similarly, and we characterize all prime spectral tetris frames. This character- 
ization ultimately allows us to determine precisely when the spectral tetris construction works 
for M/N < 2 (Corollary 021) . 

The desirable properties possessed by $2 is not only a low redundancy and avoiding repe- 
tition of frame vectors, but also that it is not the union of two tight frames as 4> is. In fact, we 
would like to argue that tight frames that are unions of tight frames (as $ above) are undesirable 
for encoding for the same reason tight frames being union of orthonormal bases are undesirable: 
they lead to encoding schemes where each sub-frame carry the same information. In particular, 
if <I> = $1 U 3»2 5 then f = g + h represented in terms of corresponds to representing g in <3>i 
and h in 4>2 completely independent and then adding the two representations in the final stage. 
This can be desirable, e.g., for image separation using sparsity techniques Il22l . Here one would 
choose 4>i and 4>2 to be incoherent tight frames so that, loosely speaking, the frame vectors of 
4>i look like noise when represented in terms of 4>2 and vice versa. However, if 4>i and 4>2 are 
very coherent, this setup is often not desirable, e.g., in the spectral tetris example the two frames 
4>i and 4>2 have maximal coherence hence taking the union of these frames will not bring us 
anything new. 

Our paper is organized as follows. In Section |2] we review some of the basic properties of 
frames. In Section [3] we define the notion of prime tight frames and analyze their structures, 
and, in particular, we prove the decomposition result, Theorem 13.41 mentioned above. In sec- 
tion HI we then restrict to prime unit-norm tight frames, and characterize the prime tight frames 
among certain families of FUNTFs including equiangular tight frames, harmonic tight frame, 
and spectral tetris frames. 

2 Notation and basic properties of frames 

In this section we collect some of the basic definitions and properties of frames. We let IK denote 
either M or C, and work in K N , N > 1, considered as a vector space over the scalar field IK. Of 
course, all results of this paper generalize from to the setting of general iV-dimensional 
Hilbert spaces. 

Definition 2.1. Given M >N,a collection of vectors 4> = {(pj}^ C K N is called & finite frame 
for ~K N if there are two constants < A < B such that 



A||x|| 2 < ^Kjc,^)! 2 ^^!^!! 2 , forallxeK^. (2.1) 



M 

E 

;=1 



If the frame bounds A and B are equal, the frame {<p,}™ =1 C IK is called a finite tight frame for 
K N . In this case, 



M 

A\\x\\ 2 = £ I (x, (ft) | 2 , for all x G K N . (2.2) 

;=1 

A finite tight frame <3> = {(p,}^ C K N consisting of equal norm vectors is called a isometric 
tight frame for ~K N . In this case, the frame bound is A = Mc/N, where c = ||<Pfc|| for all k = 
1,2...,M. 

When c = 1 and <I> = {(pj}^ C K N is an isometric tight frame, we refer to $ as afinite unit 
norm tight frame (FUNTF) for K N . In this case, the frame bound is A = M/N. 
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We shall also write 4> = [(pi]fL { £ Mat(A/ x M, K) asaniVxM matrix. Hence, we sometimes 
consider <J> as a collection of (column) vectors and sometimes as a matrix with coordinate 
vectors written as columns; this should not lead to ambiguity. We further use the matrix notation 

<t> = [<Pij] or [&]ij = qHj. 

Given a collection of M vectors $ = {(pi}^ in K N , the synthesis operator is the mapping 

M 

F : K M — > K N , (c l -)£i^*(c/) = £c i ^. 

j=l 

Its adjoint operator is called the analysis operator and given by 

F*:K N ^K M , x^®*x=((x,(pi))f =v 

Using these operators, it is easy to see that {(p,}^ is a frame if and only if the frame 
operator defined by 

M 

S := FF* : K N K N , x m- <t>®*x = £ (x, (pi) <p ; 

(=1 

is positive, self-adjoint, and invertible. In this case, the following reconstruction formula holds 

M 

x = £ (x t 5- 1 cp t ) (p h for all x G K N , (2.3) 

j=l 

and {S tyi}^ is a frame too, since S is invertible; it is called the canonical dual frame. A 
set of equal-norm vectors {<pj}^j is an isometric tight frame if and only if its frame operator 
S = AI, where A = Mc/N, c = \\(pk\\,k= 1,2, ... ,M. We mention that the frame constant A is 
referred to as the redundancy of <&. 

The Grammian or Gram operator of set of vectors 4> = {<Pf}^] C defined by 

M 

G := F*F : K M -)• K M , x ^ 4>*4>x = ( £ (<ft, (?j)xj)f^. 

7=1 

Tight frames in general, and FUNTFs in particular, have attracted a lot of attention in recent 
years due to their numerous applications. For instance, the existence and characterization of 
FUNTFs was settled by Benedetto and Fickus in [0, through the frame potential defined by 

M M 

fp(MEi) = LLKm^ (2.4) 

1=1 7=1 

In particular, they characterize FUNTFs as minimizers of the frame potential: 

Theorem 2.2 ([|2l Theorem 7.1]). Let M be fixed and consider the minimization of the frame 
potential among all collections ofM points on the sphere S N ~ [ . 

(i) IfM<N, then the minimum of the frame potential is M. The minimizers are exactly the 
orthonormal systems for ~K. N with M elements. 



4of|22] 



J. Lemvig, C. Miller, K. A. Okoudjou 



Prime tight frames 



(ii) IfM > N, then the minimum of the frame potential is The minimizers are exactly the 
FUNTFsfor K N with M elements. 

Casazza, Fickus, Kovacevic, Leon and Tremain [|9| extended this result and obtained a char- 
acterization of all tight frames as minimizers of a corresponding frame potential, under certain 
conditions on the norms of the frame vectors. The frame potential has been further generalized 
to characterize all tight frames that have an additional convolutional structure [[141 [T9l . These 
type of frames, referred to as convolutional frames, are essentially models for filter banks used 
in signal processing. We refer to lfl2ll for related results on tight frames. 

3 Definition of prime and divisible tight frames 

We define prime and divisible tight frames and prove certain of their elementary properties in 
Section 13.11 In Section 13.21 we prove that prime frames exist in every dimension. This allows 
us to prove that prime tight frames actually form a dense subset of the set of all tight frames 
(Section 1331) . Consequently, every tight frames can be approximated with a prime one. 

3.1 Definitions and elementary properties 

We now formally introduce the notion of prime and divisible tight frames. 

Definition 3.1. Let M > N be given. A tight frame 4> = {(pi\f =l in K N is said to be prime, if 
no proper subset of 4> is a tight frame for K N . If 4> is not prime, we say that it is divisible. In 
particular, given an integer p with N < p < M —N, the tight frame <E> is (M, p) -divisible if there 
exists a subset of <E> containing p vectors that is also a tight frame. 

The restriction N < p < M — N for p E N in the definition of (M,p) -divisible frames is 
motivated by the following trivial observation. 

Remark 1. For a tight frame $ = {(pi]^L x in K N to be prime it is sufficient that M < IN. In 
other words, for <I> to be divisible it is necessary that M > 2N. This follows from the fact that, 
if M < IN, it is impossible to partition $ into two spanning sets. 

We note that there exist frames <E> that are (M,p) -divisible for all p in the full range N < 
p < M — N, e.g., choosing $ to be the harmonic tight frame in C 2 with 24 frame vectors. On 
the other hand, it is not possible for a tight frame to be robust with respect to tightness for any 
p erasures where p E N is fixed, hence within the class of tight frames the erasure problem has 
no solution. This negative result and a simple symmetry observation are stated in the following 
proposition. 

Proposition 3.2. Let M >N >2be given. Let & = {(pj}f =l be a tight frame in K N . 

(i) 4> is (M,p)-divisible if and only z/4> is {M,N — p) -divisible. 

(ii) Given N < p < M — N, there exists no (M, p) -divisible tight frame 4> in K N for which 
every sub-collection of<t> consisting of p vectors is a tight frame. 

Proof, (i): The proof of this part is trivial and so we omit it. 

(ii): First assume that 4> = {(pk}^ = i is a FUNTF that is (M,p) -divisible and such that every 
sub-collection of p vectors is again a FUNTF. Let 4>i = {<p,}; e y, where J = {1, . . .,/?}, be the 
tight frame of the first p vectors from 4>. We replace the £th frame vector in 4>i by the £th vector 
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in <£, where p + 1 <k<M, and denote this tight frame $2 := {<Pi}ieJ\{£} u Then, given 
any x G K N , we have 

SlWl 2 = El^w>l 2 = E i(x,«p ; )i 2 +i(x,(p,)i 2 . 

This implies that | (jc, <j^) | = | (x, (j^) | • Since £ and are arbitrary, we have that | (x, <j^) | = | (x, 9^) | 
for each x G K N , each I — 1,2, ... ,77, and each p + 1 < k < M. In particular, we have that 
I (<Pk, <Pf)| = 1 for 1 < £ < p and p + 1 <k<M. Since the vectors <p, are unit-norm, for a fixed 
k = ko <E {p + 1,. . . ,M}, this implies that <jfy G span (^ ) for all 1 < £ < p. Hence, the span of 
4>i is one-dimensional which contradicts the assumption that $1 is a tight frame. We conclude 
that there cannot exist a FUNTF that is (M, p) -divisible and for which all sub-collections of p 
vectors are also FUNTF. This proof carries over also to the case of equal norm tight frames. 

The last case to check is when <J> = {(p,}^ is an (M,p) -divisible non-equal norm, tight 
frame such that any subset of p vectors is again a tight frame. Assume without loss of generality 
that ||<pi|| > ||<p^|| > .. . > ||<Pm||- Let 4>i = {<p,}f =1 be the subset of the first p vectors from 4>, 
and let $2 = {tyiYiZi u {<Pk} where p + 1 < k < M is the first index such that ||<p p || > \\<Pk\\- 
Then, given any x G K N , we have 

^w 2 =fi(^.)i 2 

;=1 



and 



^W N| 2 = P £ |(jC)<Pi . )|2 + |(X)(pfe)| : 



This implies that 
for each x G K N . If we choose x = <p p , we find that 



llrf-l(^^)l 2 =^™ll^ll 2 , 

and similarly, for x = (pk, we have 

I(^%)| 2 -Ilrf = ^^llrf- 
Combining these two equations we obtain 

which leads to N = 1 . This contradiction concludes the proof. □ 

The following result shows that if we can take out a subset of vectors from a tight frame 
such that these vectors form a tight frame, then what is left is automatically also a tight frame. 
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Lemma 3.3. Let M > 2N. Suppose 4> = {<Pk}f = i is a divisible tight frame for K , and let 
$1 C4> denote a tight frame for ¥± N . Then 4>2 = \ 3>i is also a tight frame for K. N . 

Proof. Assume that and $1 C$ are tight frames with frame bounds A and A\, respectively. 
Note that A\ < A. Let <J>2 = 4>\4>i. After a possible reordering of columns, we have that 
4> = [4>i $2]. Since 4> and 4>i are tight frames, it follows that 



*2. 



<D 1 $*+$ 2 <D*=A 1 / JV + $ 2 ^ 1 (3.1) 



AI N = = [4>i <J> 2 ] 
which implies that 

<$> 2 <$>* 2 = {A-A l )I N . (3.2) 
Hence, we see that 4>2 is tight frame with frame bound A — A\ > 0. □ 

Lemma 13.31 is a "reverse" to the well-known fact that the union of tight frames is a tight 
frame. We say that tight frames 4>i and 4>2 as in Lemma 1331 are divisors or factors of 4>. We 
refer to Definition [33] below for more details. 

The main observation of this section is that every tight frame of M vectors in K N is a finite 
union of prime tight frames. 

Theorem 3.4. Given M > N, every tight frame of M vectors in K is a finite union of prime 
tight frames. 

Proof. If 4> is prime, there is nothing to prove. So suppose that is divisible tight frame with 
frame bound A. Then, by definition, we can partition 4> into 4> = 4>i U (4> \ 4>i ) , where $i C$ 
is a tight frame. It follows that 4>2 := 4>\4>i is also a tight frame by Lemma [331 If 4>i and 
4>2 are prime, we are done. If not, repeat the process. In each step of this procedure we split 
a tight frame into two sets of cardinality at least each. Hence, by Remark [U the procedure 
terminates after at most [JCJ steps. □ 

The decomposition result in Theorem 13 .41 naturally leads to the following definition. 

Definition 3.5. Let 4> be a tight frame. Suppose, for some K e N, 

4> = 4>!U---U^, (3.3) 

where 4> fc , k = 1 , . . . ,K, are prime tight frames. We shall say that <i> k are prime factors or prime 
divisors of 4>. 

The factors are in general not unique as the following examples show. 

Example 1. (a) The 6th roots of unity FUNTF in M 2 , that is, the frame of six unit-norm vectors 
each 2n/6 apart as the vertices of a hexagon, factors into two FUNTFs, each of which 
consists of three vectors. But these are not unique since you can choose the two set of 
three frame vectors in eight different ways. However, up to multiplication of individual 
frame vectors by —1 this decomposition is in fact unique. Moreover, you cannot factor this 
FUNTF into prime FUNTFs of, say, size not equal to three. In particular, you cannot factor 
it into three FUNTFs of size two. 
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(b) The uniqueness of the factors up to scalar multiplication in part (a) of the example does not 
hold in general. The harmonic frame of 10 vectors in C 2 decomposes into either five frames 
of size two or two frames of size five. We refer to Section for more details on this. 

Tight frames 4> and *F are unitarily equivalent if there is a bijection p : {1,...,M} — > 
{I,.. . ,M}, a unitary N x N matrix U and ae,-GK with |c,-| = c > such that Xj/f = ciU (p p ij\ for 
all i = l,...,M. In matrix notation this reads *P = U<i>PC, where P is the M x M permutation 
matrix for p and C = diag (c\ , . . . , cm)- Note that there exists other related equivalence notions, 
e.g., in IfTOl c, is replaced by a fixed positive scalar c > 0. We remark that it is only necessary 
to introduce permutations in the equivalence relations if one considers a frame as a sequence. If 
a frame is considered as a non-ordered collection of vectors (with repetition allowed), then this 
is clearly not necessary. We refer to lfT8ll for more on equivalence relations between FUNTFs. 
The following result shows that prime frame are equivalence classes. 

Proposition 3.6. Suppose <I> and *P are unitarily equivalent tight frames for K N . Then $ is 
prime if and only if*¥ is prime. 

Proof. We prove that <J» is divisible if and only if *F is divisible. Assume 4>y := [<?>,•];<=./ is a 
divisor of 4> for some J C {1, . . . ,M}. Then Qj®} = AI N for some A > 0. Since Vp-Ui) — 
c p -iffdJ (pi by assumption, we have that 

where Cj = diag(Cp-i^) . This shows that *P is divisible with divisor M^-i^. By the sym- 
metry of the equivalence relation, this completes the proof. □ 

Proposition 13.61 also holds for the notion of unitarily equivalence used in IflOl . However, the 
result is false if one uses the coarser notion of equivalence where the matrix U is only assumed 
to be invertible. 

Remark 2. The notion of prime tight frames easily generalizes to non-tight frames. Recall that 
if 4> is a frame with frame operator S = and Grammian G = $*4>, then the associated 
canonical Parseval frame is given by 5 _1//2 4>. We shall then say that $ is a prime frame if its 
canonical Parseval frame is a prime tight frame. One can describe such prime non-tight frames 
4> in terms of certain properties of eigenbases of the Grammian. Let {v;} ;=1 be any orthonormal 
eigenbasis associated with the non-zero eigenvalues of the Grammian G. Set V* = [vi]f =l & 
K MxN . Then 4> is prime if and only if the tight frame V is prime. In other words, a frame 4> 
is prime precisely when one can find a orthonormal eigenbasis {v,} ;=1 C K. M associated with 
the non-zero eigenvalues of the Grammian G that remains orthonormal (up to a scaling factor) 
after projecting it onto a coordinate hyperplane of K M . To see this claim, let 4> = UYW* be 
the full singular value decomposition of where U is an N x N unitary and W is an M x M 
unitary. Further, let $ = ULW* be the compact form of the singular value decomposition, where 
E = diag (oi, . . . , On) are the singular values of and W the first N columns of W. A direct 
calculation shows that 5 ,_1 / 2 <J> = UW*, hence, by Proposition 13 .61 the canonical tight frame is 
prime if and only if W* is prime. Since 

G = = WL*U*ULW* = WAW* = WAW*, 
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where A = diag (of, . . . , o N ) and A = diag (of , . . . , 0$, . . . , 0) G K MxM , we see that {vi}f =1 
defined by V* = [vj]f =l = W is an eigenbasis associated with the N non-zero eigenvalues of the 
Grammian G. This proves the claim. 

If all non-zero eigenvalues are distinct, V is uniquely determined up to multiplication by 
scalar with modulus one. If some eigenvalues has algebraic multiplicity greater than one, then 
V is no longer uniquely determined. However, if one such V is prime so are also all other 
choices of V since they are unitarily equivalent. 

3.2 Existence of prime tight frames 

Let 4> be a union of an orthonormal basis for K N and M — N zero vectors. Then <3> is prime 
for any N,M E N. This trivial observation shows the existence of prime tight frames for all 
M > N E N. The following result extends this fact to tight frames of non-zero vectors. 

Proposition 3.7. For each dimension N EN there exists a prime tight frame for ¥L N with M 
non-zero vectors for any M > N 6 N. 

Proof. The case N — 1 is trivial, hence we assume N >2. Let *P E K^- 1 )^- 1 ) be a Parseval 
frame for K^ -1 . We now extend *P into an N x M matrix by first adding a 1 x (M-l) row 
vector with zeros as a new Mh row, and then by adding e^ E ~K N as a new Mth column. The 
picture is: 

'•• : ••' 
... vj/ ... o 

c£> = 

••' : '■• 
1_ 

The new frame 4> is a Parseval frame for K N since 4>4>* = In- Moreover, it is prime since any 
tight frame 4>i C$ must contain the Mth column in order to span ~K N . □ 

Remark 3. For certain values of N,M, we can extend Proposition l3.7l to show existence of prime 
unit-norm tight frames. The case N < M < IN follows from Remark [T] If M > 2N and M is 
prime, then we take the M x M DFT matrix and choose any Af rows to be our frame <I>. Since 
no proper subset of the (primitive) Mth root of unity sum to zero, there exists no way to divide 
4> into two FUNTFs. Hence, $ is a prime FUNTF. 

The existence of prime tight frames allows us to prove that most tight frame are prime. 
Precisely what is meant by "most" is the topic of Section 1331 We start here by showing that a 
random tight frame is prime with probability one. 

Let V be an N x M matrix with entries independently drawn at random from an absolutely 
continuous distribution with respect to the Lebesgue measure; a standard choice could be the 
Gaussian distribution of zero mean and unit variance. With probability one, V is a frame and 
thus S- l / 2 V is a Parseval frame, where S is the frame operator of V. Alternatively, we can 
obtain a Parseval frame from V by performing the Gram-Schmidt algorithm on the rows of V 
or by keeping iV rows of a random orthogonal M x M matrix. We call such frames for random, 
tight frames. As one would expect random tight frames are prime. 

Theorem 3.8. Let N,M EN. Random tight frames of M vectors in M. N are prime with proba- 
bility one. 
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Proof. Without loss of generality we assume the random tight = {(pi} i=1 to be constructed 
as follows. Let V be an N x M matrix with entries independently drawn at random from a 
absolutely continuous distribution. We denote the rows of V by Vj, j = 1 , . . . , N. Let the rows Uj 
of the matrix U be given as follows: uj = Vj — Y,{Zi P r oj« A . (v/)» where the projection of v onto u 
is defined by proj ;( (v) = (v,u) \\u\\ u. Finally, we normalize the row vectors uj, j = 1,. .. ,N, 
to have norm A 2 for some A > and let these normalized vectors be the rows of 4>. Then <3> is 
a tight frame with bound A. 

Assume towards a contradiction that 4> = {(pj}f =l is divisible, and let 4> = {<P/}, e ; and 
<J> = {(Pi} ie j be divisor of 4> = 4> U 4>. Let v 7 and Vj be the row vector obtained from vj by 
keeping the indices in / and /, respectively, and similarly for the row vector Uj. To obtain a 
contradiction, it suffices to show that ii\ and w 2 are not orthogonal. Since u\ = vi, we see that 

(V2,Ml) _ \ ,~ -\ (V1,V2) /- -v 
(Ml, "2) = (Vl,«2) = ( Vl,V2-— -2~Wl ) = (Vl,V 2 )-- -2-(Vi,Vi) 

l M l|| / ll v l|l 

- i|2 n~ ii2 

Vl,V 2 )-- -2 (Vi,V 2 ). 



1 \\1 s LJ z/ 11 ll2 
l v l|| Vi 



Hence, u\ and m 2 are orthogonal if and only if 

(vi,v 2 ) (vi,v 2 ) 



l~ l|2 ' ||" ||2 ' 
l v l|| vi 



By Proposition 13 .71 not all random vectors vi and v 2 will satisfy this equation. Therefore, since 
vi and v 2 are random vectors, this algebraic equation (written in real and imaginary parts) is 
with probability one not satisfied, see Section [331 □ 



3.3 Prime tight frames are dense 

By the fact that the union of any two tight frames again is a tight frame (and, e.g., by Propo- 
sition [X6]) we know that there exists many divisible frames. However, we shall show that in 
a strong sense there are very few divisible frames within the set of A-tight frames, where A is 
fixed. To state this result, we need to set some notations. 

For A > fixed, let ^(N,M,A) be the set of all tight frames with frame bound A, that is, 

&{N,M,A) = e Mat(N x M, IK) : = AI N } . 

Let 

&>(N,M,A) = {<£ e &(N,M,A) : 4> is prime} . 

When A = 1, we simply denote these sets JP(N,M) and ^(N,M). Note that JP(N,M) is the 
Stiefel manifold, see [281, which is invariant under multiplication by U E U (N) from the left and 
by U E U (M) from the right. Consequently, there exists an invariant Haar probability measure 
/l on ^(N,M). The results in this section hold for any measure that is absolutely continuous 
with respect to /i. 

Following the setup in [1 , Section 3], each of the entries in a matrix 4> E J^(N,M) is written 
in terms of its real and imaginary parts: <p ; j = Xij + iytj for all 1, j. This way we can consider 
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JP(N,M) as an real algebraic variety in M . Since ^(N,M) is an irreducible variety, every 
non-empty Zariski-open subset of JP(N,M) is open and dense in ^(N,M) in the (induced) 
standard topology 0]|. Moreover, the complement of a non-empty Zariski-open subset is of 
[i -measure zero. The following result says that being prime within the set of tight frames is a 
very generic property. 

Theorem 3.9. LetN,M e N and A > 0. The set of prime A-tight frames &(N,M,A) is open and 
dense in the set of all A-tight frames ,M,A) in the (induced) standard topology. Moreover, 
the complement c?(N,M,A) \ &{N,M,A) is of measure zero in 11. 

Proof. We can without loss of generality take A = 1. Let <I> £ ^(N,M) and let S be the power 
set of {1, . . . ,M}. We use the notation = [(pi\iei for I E S. Now, $ is divisible if and only if 

= cIn 

for some ^ / £ S and c > 0. These orthogonality conditions can be expressed as polynomial 
equations in xij and yy introduced above. Since S is finite, the set of prime 1 -tight frames 
&(N,M) is a finite intersection of the complement of such sets in J^(N,M). The sub-variety 
&(N,M) is therefore Zariski-open in the irreducible variety ^(N,M). By Proposition 13.71 the 
set &(N,M) is non-empty. Since &(N,M) is a non-empty Zariski-open set in an irreducible 
variety, the result follows. □ 

We remark that being Zariski-open is a strong property even in a measure theoretic sense: 
being open and dense in the standard topology does not imply having a complement of measure 
zero, but being Zariski-open does imply having a complement of measure zero. 

By Theorem 13.91 we see that if 4> is a (divisible or prime) tight frame and <t> a random, 
arbitrarily small perturbation of <I> such that <t> again is tight, then $ is prime with probability 
one. From Theorem 13 .91 we also have the following density result. 

Corollary 3.10. Every tight frame is arbitrarily close to a prime tight frame. 

4 Classification of certain prime tight frames 

In this section we characterize prime frames within well-known families of FUNTFs. More 
specifically, we consider equiangular tight frames, harmonic frames, and spectral tetris frames. 
Since generalized harmonic FUNTFs are unitarily equivalent to harmonic FUNTFs, we also 
obtain a classification for this class of FUNTFs . We refer to (El for more on generalized 
harmonic FUNTFs. 

4.1 Equiangular FUNTFs 

Equiangular tight frames are FUNTFs for which the modulus of the inner products of each pair 
of distinct frame vectors has the same value H29U30H . We show that, when they exist, equiangular 
tight frames are automatically prime tight frames. This is a necessary condition which to the 
best of our knowledge is new. Equiangular tight frames are important in many applications in 
signal processing, and they arise in many contexts such as spherical designs lfT3ll26ll . quantum 
computing Il25ll . and equiangular lines [7]. Moreover, we mention that equiangular tight frames 
are exactly sequences that achieve the Welch bounds with equality ll3"Tll3~2Tl . 
We first recall the definition of equiangular tight frames. 
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Definition 4.1. Let 4> be a FUNTF for K N . The frame = {<p,} /=1 is said to be equiangular if 
\(<Pi,<Pj) \ = c for all i,j= 1, . . . ,M with i ^ j for some constant c > 0. 

The constant c cannot be arbitrary and in fact is related to N, and M is a precise manner. 
This follows from the Welch bounds ll3Tll3~2ll : 

Theorem 4.2 ( 11291 Theorem 2.3]). Suppose $ = {(Pi\f = \ be a unit norm frame in K N . Then 



max 1(^,^)1 > 



M-N 

N(M- 1)' 



and equality holds if and only if 4> is an equiangular tight frame. Moreover, the following 
assertions hold: 

(i) ifK. = C, equality can hold only ifN < M < N 2 , 

(ii) ifK 



t, equality can hold only ifN <M< Mj^hli 



The existence of equiangular tight frames is still an open problem. However, we show below 
that, when they exist, equiangular tight frames are also prime. 



Theorem 4.3. LetN > 2. Equiangular FUNTF s ofM vectors in K , when they exist, are prime. 

Proof. Assume that 4> = {(ft}^ C K is an equiangular tight frame. Assume towards a con- 
tradiction that 4> is (M, p) -divisible for some N < p < M — N. Then write $ = $iU<I>2 where 
4>i and 4>2 are divisors of 4> of size p and M — p, respectively. Consequently, we see that <I>, 
4>i, and 4>2 are all equiangular tight frames. According to Theorem 14 .2[ 



M-N 
N{M-l) 



p-N 
N(p-l) 



M-p-N 
N{M-p-l) 



with N < p < M — N. But a series of easy calculations leads to a contradiction. Thus, 4> cannot 
be (M,p) -divisible for any N < p < M -N. □ 



The trivial case A/ = 1 is excluded from Theorem 14.31 for a good reason: equiangular 
FUNTFs in K 1 are unions of scalars of modulus one. These are, of course, always divisible. 

Finally, we remark that Grassmanian ( O) tight frames are not, in general, prime tight 
frames, e.g., any Grassmanian frame of four frame vectors in IR 2 is (4, 2)-divisible. 

4.2 Harmonic frames 

In this section we completely characterize all harmonic FUNTFs that are prime, and for those 
that are divisible we describe their factors, but first we recall the definition of harmonic frames. 
Given M > N, a harmonic tight frame (HTF) is obtained by keeping and renormalizing the 
first coordinates from anMxM discrete Fourier transform, that is, <t> = {(pk\k=o ^ s a HTF, 
where the (k+ l)th column of <I> is given by 



/ 1 \ 



(pk 



WnJ 



gC 



N 
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where s > 0, Jm '■= exp(2ni/M) is the Mth root of unity and (£>„ := y^ 1 = e 2m ( n -i)/ M f or 
n = 1, . . . ,N. When there is no risk of confusion we shall simply write y for Jm = ex P (2ni/M). 
We denote the obtained HTF by HTF(Af,M, s), and we see that this tight frame has frame bound 
A = sM/N and frame vector norms ||<p&|| = s. Hence, for s = 1 we have a unit-norm, M/N- 
tight frame, while we for s — N/M have a Parseval frame. When nothing else is mentioned we 
assume s = 1 for simplicity. 

Let us fix some notation and assumptions for this section. We will always assume N > 2 
since, if N = 1, any HTF is (M,p) -divisible for all p. Throughout this section we will denote 
the index set {1,2,..., M} by /, and we will let J\ denote a subset of / and put fa '■ = I \ Ji ■ If d 
is a divisor of M, we define index sets I(d,q) C / = {1, 2, . . .M} as follows 

l(d,q) = yc^- + q:k = 0,l,...,d-l^, q=l,...,M/d. (4.1) 

These index sets are a disjoint partition of / = {1,2,... ,M} for a fixed divisor d. For any n = 
1,2, ... : N— 1, the set {(/*)'" : m £ I{d, 1)} is a subgroup of {(/")'" : m £ /} in the circle group, 
and {(Y 1 )" 1 '■ m e K^^)} is a coset. Furthermore, we assume the following prime factorization 
M — p^p® 2 ■ ■ ■ Pr r with a i ^ N and p,- prime and ordered such that p\ > pi+\. 

Let = {<P/t}£li cC^bea HTF with index set / = {1, 2, . . . , Af}. Since $ = {(j^W is a 
tight frame, the rows of 4> are equal-norm and orthogonal. In particular, we have 

M 

0= (<p'\ <?"') = £ (y(»-i)-(»'-i)yn-i = j^yi-rfyn-i f rn^n'£{l,2,...,iV} (4.2) 

m=l m£/ 

where <p 7 denotes the j'th row of Now, $ is divisible exactly when it is possible to split the 
sum over m £ / into two sums, each summing to zero, for each n^n 1 . Here we have used that 
the norms of the Af rows of any sub-collection of $ are automatically equal since the entries 
of $ are equal in modulus. Therefore it is only necessary to consider the row-orthogonality 
requirement of potential divisors. We will use the following result repeatedly. 

Lemma 4.4. Let N,M > 2 be given, and let = {(pk}f =l C C N be a HTF. Then 4> is (M,p)- 
divisible if and only if there exists J\ C / = {1, . . . ,M} with #J\ = p such that 

0= E^ 1 )'^ 1 forne{2,...,N}. (4.3) 

me7i 

Proof. Let $1 := <!>/, = {(Pk}keJi> anc * ^ et W ^ e me ^ tn row °^ < ^ >1 ' ^ ote tnat tne e£ l uat i ons 
in (14.31) are equivalent to the statement that {q>\ , <p") = for n = 2, . . . ,N. Assume (14.31) holds. 
Then, since y~ k = y M ~ , we see that 

0= £ (y"-«')'«-l 
me/i 

holds for all n, n' £ {2, . . . , A 7 } with n 7^ n ; . The last statement is equivalent to (<p", <p" ) = for 
all n^«'e{2,.. . , N} which in matrix notation becomes <i>i<i>\ = cl n . Therefore, equation (14.31) 
implies that $1 is a tight frame. The converse implication follows easily from the above. □ 
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Suppose is (M,p) -divisible. By taking n = 2 in (14.31) we see that the sub-sum over Ji 
must be the sum of p Mth roots of unity, and, of course, the second sum over J2 must be the 
sum of (M — p) Mth roots of unity. This is an example of vanishing sums of roots of unity Il24l . 
When one of the vanishing sub-sums contains p distinct Mth roots of unity, one says that M is 
p-balancing ll27l . In particular, we shall use the following result proved in lETl . 

Theorem 4.5 ( fl2Vl Theorem 2]). Let M = p^p" 2 ■■■p? r with each Pi prime and each a, > 0. 
Then M is k-balancing if and only if both k and M — k are in Nopi + N0P2 + • • • + ^oPr, that is, 
both k and M — k are linear combination of p\ with nonnegative integer coefficients. 

By Lemma I4~4l we immediately have that, for N = 2, a HTF is (M,p) -divisible if and only 
if M is p-balancing. Therefore, Theorem I4.5l is enough to characterize all prime HTFs in C 2 . 

Theorem 4.6. Let M > 2 be given. Suppose 4> = {(Pk}f =l C C 2 is a HTF. Then 4> is prime if 
and only ifM is a prime integer. Furthermore, if M is not prime, then 4> is (M ', d) -divisible for 
each 2 < d < M — 2 for which both d andM — d are in f%pi + No/>2 + ■ • • + ^oPr, in particular, 
$ is (M, d) -divisible for every divisor d ofM. 

Proof. We will use the fact that the frame HTF(2,M,s) is prime if and only if M is not d- 
balancing for any 2 < d < M — 2. 

Assume first that $ is prime, that is, that M is not J-balancing for any 2 < d < M — 2. 
Towards a contradiction assume further that M is not prime so that M = p" 1 1 p^ 2 ■ ■ ■ p® r with 
r > 1. Then M = p\b with b — pf 1- p% 2 ■ ■ ■ p" r > 2. For d = p\ we have 

M-d = pib-p l =p 1 (b-l), 

but this contradicts the fact that M is not J-balancing for any 2 < d < M — 2. 

For the opposite implication, we observe that the prime factorization of M = p\ is trivial 
when M itself is prime. Since there is no divisor < d < M so that d,M — d E MNq, we see that 
M is not J-balancing for any N < d < M — N. Thus, we have proved the bi-implication part of 
the theorem. 

The "furthermore" statement follows immediately from the above and Theorem 14 .51 □ 

Remark 4. The bi-implication part of Theorem l4.6l can also be proved without the use of Theo- 
rem !4.5l as follows: 

If M is prime, then = e lmk / M is a primitive Mth root of unity for each k= 1,2, . . . ,M— 1. 
Hence, no proper subset of these roots of unity can sum to zero. Consequently, 4> is prime. On 
the other hand, if 4> is prime, and if M has a divisor d > 2, then M = bd for some b > 1. If one 
takes J 1 = {nb : n = 1, . . . , d} = {b,2b, .. .,db = M}, then <t>j t is a tight frame by Lemma l4~4l 
which is a contradiction. 

In general, when N > 2, the characterization of prime HTFs is a little more involved. Indeed, 
we now have multiple rows consisting of Mth roots of unity, and we must insure, for (M,d)- 
divisibility, not only that each of these row is J-balancing, but also that the subset of d roots 
that sum to zero in each of these rows, comes from the same columns. 

In order to formally state this we need to define a few sets of integers. Using the notation 
fixed above, we define 

D m ,n = {d G {N,...,M -N} : d is a divisor of M}, (4.4) 
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Pm,n = Dm,n \ {d G Dm,n ■ 3c G Dm,n such that c is a divisor of d} , (4.5) 

and 

K K 

Sm,n = {se {N,...,M-N} :s = Y,a k qk,M-s = ^hqk, where a kl b k G N 0l q k G Pm,n}, 

k=l k=l 

where A" = #Pm,n- Note that 5 G Smjv if and only if M — 5 G Sm,n- It is also clear that 

and that Dm^ is empty for any N G N if M is prime or if M < 2N. Note that if M > 2N, the 
condition d < M — N in the definition of Dm,n is redundant as no divisor of M can be greater than 
M — N. Moreover, the set Sm,n is empty, if and only if Dm.n is empty. The significance of these 
sets will become evident in Theorem 14.81 and 14.101 and Corollary 14.91 below, but we mention 
here that Pm,n will determine the cardinality of the prime factors and Sm.n the cardinality of 
every possible divisor of the HTF. In the following example we calculate these sets for various 
M,NeN. 

Example 2. (a) For prime M G Z and any N G N, we have Pm,n = L>m,n = Sm,n = 0- 

(b) For M = 9, TV = 2 or TV = 3: D M , N = Pmjj = {3}, S m ,n = {3,6}. 

(c) For M = 9,N > 4: D M , N = P m ,n = S m ,n = ©• 

(d) For M = 10, N = 2: D M:N = P m ,n = {2,5}, S m ,n = {2, 4, 5, 6, 8}. 

(e) For M = 10, /V = 3,4 or 5: D m ,n = Pm,n = S m ,n = {5}. 

(f) For M = 24, N = 2: D m ,n = {2, 3, 4, 6, 8, 12}, P m ,n = {2, 3}, S m ,n = {2,3,..., 22}. 

(g) For M = 24, N = 3: D M , N = {3, 4, 6, 8, 12}, P M ^ = {3, 4}, S m ,jv = {3, 4} U {6, 7, . . . , 1 8} U 
{20,21}. 

(h) For M = 24, /V = 4: D M ,# = {4, 6, 8, 12}, P M ,w = {4, 6}, S M jv = {4,6,..., 20}. 
We will also need a standard result from Hewitt and Ross ifTTl . 

Lemma 4.7 ( lfT7l Lemma 23. 19]). Let d G N be given. Then, for m G Z, 

1 ^ » .few f 1 m c ,7Z 

I Y e 2ni ^f = 1 ' (4.6) 

d^ [0 meZ\dZ. 

We are now ready to state and prove the following result is a characterization of HTFs that 
are prime. 

Theorem 4.8. Let M = p^p® 2 • ■■pf r with pi prime and a,- > for i = 1, . . . ,r. Assume that 
= {'Pa:}^! C is a HTF. Then $ is prime if and only if Dm,n = 0- In particular, if M is 
prime, then the frame 4> is prime. 

Proof. Assume that 4> is prime and let us prove that Dm,n = ®- Assume by way of contradiction 
that there exists d G Dm.n 7^ 0, i.e., J > iV is a divisor of M. Take 7i = /(J, 1), where I(d, 1) is 
defined in (14.11) . Then, for n G {2, . . . , J}, 

I (iSr 1 )™- 1 = I eTTfc-WT) = £ e 2m ^ = 0, 

m€/i ifc=0 fc=0 
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where the last equality follows from Lemma 14771 Since d > N, we see by Lemma 14741 that 4> is 
(M, d) -divisible with divisor 4>i := <J>/j which is a contradiction. 

We now prove that if 4> is not prime, then D m ,n ®- Thus assume that 4> is (M, J) -divisible 
for some d such that N < d < M — N. Let 4> = 4>i U 4>2 be divisors of 4>, where 4>i = ^ 
and $2 = &j 2 w i m #A = d and #72 = M — d. By Lemma 14741 we see that our assumption is 
equivalent to assuming the existence of a index set J\ C {1, . . . ,M} of cardinality d such that 

E dfi" 1 ) 1 "" 1 = I ^ 1= (4-7) 

holds for each n = 2, . . . ,N, where J" := (n — 1)J\ mod M. In particular, for n = 2, this means 
that M is -balancing, hence d = Y*l=l a kPk an d M — d = YJ k =\ bkPk for a^bk G No- Let i?i C 
{ 1 , . . . , r} be the indices for which a& 7^ 0. We see that the index set J\ has the following form 

7 i = U U I (Pk,<lk)= IJ h, where 4:= [J I(p k ,qk), 
keRiq k eQ k keRi q k eQ k 

for some Qk C {1, . . . ,M/ /)^} with #<2/t = for each k E R\. We can assume without loss of 
generality that, for 1, j £ Ri, 

KiPi 7^ f° r each rii = 1 , . . . , a,- and nj = 1 , . . . , a,-, (4.8) 

whenever i 7^ j. 

We need to prove that Dm,n 7^ 0- It follows from d = Y^keR\ a kPk that, if p^ > iV for some 
ko E R\, then pk G £>m,n 7^ 0, which concludes the proof. If however, pk < N for each k ER\, 
we claim that > N for each k G To show this claim, let &o G i?i be fixed and assume 
Pk < N for all k E R\. Since pk < N is, prime, it follows that 

for n = 2, . . .,Pk,Pk + 2, . . . for k E R\. We remark that is not a ^-multiple of any of the 
other prime numbers pj for nj = 1, . . . , ay by equation (14.81) which, in turn, implies that 

E (7m )" 7 - 1 =0 

m€Ji\k Q 

for n = 7?£ + 1 • Hence, it follows by the fact that 4>i is a tight frame, that 

E (YM°) m - l = Pk E (r5°r 1 = o- (4-9) 

The last equality is only possible if the index set Qk is of a certain size, in particular, #Qk > 2 
must be a divisor of Mj pk Q . If M/ pk is prime, then = #Qk = M / pk which is impossible. 
Hence, M/pk is not prime. Let d\ be the smallest divisor of M/pk . By (14.91 ) we then see 
that #Qk = nd\ for some n E N since must be a union of rc sets, say D], £ = 1, . . . ,n, of 
cardinality dj. If d\p^ > N, we are done. If not, we consider row n = d\pk + 1 in (14.71) . 
Repeating the argument above leads to 

E (YM Pk °) q - l =0. (4.10) 
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Note that (7^ Pk °) q ~ 1 is same Mth root for all q G D\ when £ is fixed. If Mj (dipk ) is prime, 
then, by (14.91) and (14.101) . #Qk = Mf which is impossible. Hence, Mj {d\p^) is not prime, 
and we let J2 be the smallest divisor of Mj (dip^). By (14.91 ) and (14.101) we then see that #Qk — 
nd\d2 for some n G N. We can continue the argument which proves the claim. 

By the proof of the claim, we also see that ak Pk is a divisor of M. Since we just proved 
a k Pk > N, we arrive at the conclusion a^p^ G D m ,n ®< 

The proof of the last statement in Theorem 14 .81 is an easy consequence of the fact that if M 
is prime, then D m ,n = 0- 

□ 

Remark 5. (a) The last part of the above result states that a HTF <E> = {(Pk\k=\ c ^ w ^ m 
M > IN is prime if M is prime. The converse of this is not true. Indeed, consider the HTF 
with M = 9 and N = 4, see Example Etc). This frame is prime, but M is not. 

(b) Theorem [43] classifies all k for which M is ^-balancing, but it does not give how to choose 
the k distinct roots out of the M roots of unity. What the proof of Theorem [478] shows is that 
for a HTF to be (M, d) -divisible, in addition to having a vanishing sub-sum of d distinct 
Mth roots of unity, we must also ensure that the sub-sum of every nth power of these roots 
of unity vanishes for n = 2, . . . , N — 1. We refer to Theorem 14 . 9 1 below for a statement on 
how to choose the d distinct roots of unity so that all these sub-sums vanishes. 

In case the HTF is divisible i.e., Dm,n 7^ 0, the following result tells us how to factor it into 
prime divisors. 

Corollary 4.9. Let N,M e N and s > be given, and let $ = {(p k }f =l = HTF{N,M,s) C C N . 
Suppose that Dm.n 7^ an d d < M. Then the following assertions are equivalent: 

(i) deP M ,N, 

(ii) 4> factors into M/d prime FUNTFs each of cardinality d, that is, 

M/d 

<*>= u*«- 

with <t>i being prime and #<&; = d. 

Furthermore, one of the factors from assertion (ii), say <&i, is the HTF with d frame vectors is 
C N , that is, $1 = HTF(Af, J, s). Let U = diag (y°, y 1 , . . . , y^ 1 ) E U(N), where y = exp (2m/M). 
The other factors are then given as <I> ; - = U l ~ l Q\for i = 1, . . . ,M/d. 

Proof. This follows from the proof of Theorem 14 .91 □ 

By Corollary 14.91 we know exactly how the prime "building blocks" of a divisible HTF 
look, hence we can use this structure to build larger divisors of the HTF. Suppose <I> is a HTF 
and d G Dm,n 7^ 0- Then 4> is (M, d) -divisible, and, moreover, {<Pi}i e iUq) 1S a tx S n t frame for 
any q = 1, . . . ,M/d. Now, we can combine these M/d tight frames of cardinality d into tight 
frames of cardinality d,2d,3<i, . . . ,M. Hence, 4> is also (M, nd) -divisible for n = 1, . . . ,M/d— 1. 
Assume further that M has another divisor which are also greater than N, say d G Dm,n- We 
can then combine unions of {(Pi} iel ^ 9 ) with unions of {<P/} ;e /(j^), where q = 1, . . . ,M/<i and 
g = 1,. . .,M/d, respectively, as long as the index sets are disjoint. Hence, to combine tight 
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frames from different divisors, we need to make sure that the same frame element <p,- does not 
appear in both frames. We make these observations precise in the following result. 

Theorem 4.10. Let M > N > 2, and 4> = {<h}k=\ C C N be a HTF. If D m ,n ^ ©, then 4> is 
(M ,s)-divisible for each s G Sm,n- 

Proof. If D m ,n 7^ 0, then is (M, J)-divisible for any divisor d of M such that N <d <M-N. 
We now show that 4>, in fact, is (M,s) -divisible for each s G Sm.n- By symmetry of the set 
Sm,n, it suffices to show that 4> is (M, ^-divisible for s G Sm,n with s < M/2. For s G Sm,n we 
have s = Y*k=\ a k<lk with aj.eNo,ft£ Pm,n and ^ = #Pm,n- For each = 1 , 2, . . . k construct an 
N x matrix <I>i based on the first g^th roots of unity. By shifting this matrix a^ — l times, 
we will have a collection of a^ such matrices. Next define 



Now is an A/ x a^^yt matrix whose rows are mutually orthogonal. We then obtain an N x d 
matrix 

which is a FUNFT. Hence, 4> is (M, s) -divisible. □ 

We end this subsection by remarking that the results of this section extend to generalized 
harmonic frame introduced by Casazza and Kovacevic [0. In (H it is shown that generalized 
HTFs are unitarily equivalent to HTFs hence, by Propo sition 13 . 6 1 and the results obtained in this 
section for HTFs, we obtain a classification of prime and divisible generalized HTFs. 

4.3 Computational aspects of divisible harmonic frames 

For harmonic tight frames the prime building blocks are exactly described by Corollary 14 .91 We 
wish to suggest a strategy that can be used to design efficient numerical tools for signal pro- 
cessing with divisible HTF. Recall that, in analyzing signals with any type of frame, one needs 
to compute c = {(x, <Pi)}.- =1 for signals x G C^. A naive way of computing the analysis step for 
divisible harmonic tight frames would be to zero pad x G C N into a vector x, G C M and then com- 
pute a full FFT of x of size M. However, Corollary 14.91 suggests a more effective analysis (and 
synthesis) process. Suppose <E> = {(pi}f =l is the divisible harmonic tight frame for C N obtained 
by keeping the first N rows of an M x M discrete Fourier matrix. Suppose p G Pm.n, where the 
set Pm,n consists of the cardinalities of the possible prime factors which is a small subset of the 



divisors of M (see equation (14.51) for the exact definition). For e.g., M = 24 and N = 4 we can 
take p = 4. The indices I q :=I(p,q) = {kM/p + q : k = 0, 1, ...,p — 1} for each q = \.. ..,M/ p 
are then the index sets of the M/ p prime factors 3>/ ? := {(f>i\ ieI . By Corollary 14 . 91 these prime 

factors are closely related, in fact, 4>/ 9 = U q ~ l &i x , where U = diag (y°, y 1 , . . ■, n ^^ 1 ) G U(N) 
is a unitary, diagonal matrix and y the Mth root of unity. Note that <t>i 1 is first N rows of 
a p x p discrete Fourier matrix. Therefore, the analysis process of computing {{x,(pi)} i=1 
for some signal x G C N can be performed by first computing y q := (U*) q ~ l x G C N for each 
q = 1, . . . ,M/p, then computing ci q := <£/ y q G C p forq= 1, . . . ,M/p, and finally combining 
ci for each q = 1, . . . ,M/p into c. The first step (U*) q ~ l xh fast and stable since U a is diagonal 
unitary, the second step is M/p FFTs of size only p, and the last step is simply a rearrangement 
of the coefficients. A similar speed-up procedure can be used for the synthesis process. We 
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think that these speed-up methods should be known in the literature since the main idea used is 
very simple and somewhat similar to the idea behind the FFT. However, we have not been able 
to find any references on this matter. We shall elaborate more on this aspect of prime harmonic 
frames elsewhere. In any case, this example illustrates that the prime factors are the key build- 
ing blocks of the frame and that one can do computations on each building block in a parallel 
and efficient way. 

4.4 Spectral Tetris frames 

Spectral tetris is a method to construct FUNTFs of M vectors in M. N with frame bound X =M/N 
for any N,M G N satisfying M > 2N. The method was developed in [0 and extended in Q. 
Here, we shall consider spectral tetris frames constructed from the algorithms in flMH under the 
usual sparsity setup that the "tetris blocks" are of size lxl and 2x2. 

We define K = {k„ : n = 0, 1 , . . . , gcd {N,M}}, where k n = ^TaW For N > M eN such 
that X = M/N > 2 the spectral tetris frame STF(iV,M) in M N is then given as the M columns of 



1 ••• 1 a\ a\ 

bi -b x 1 



mi times 



where a 



j ■' 



2 and bj 



1 a.2 Cl2 
b 2 -b 2 1 



ni2 times 



••• 1 a n -\ a n -\ 

b n -i -b n -\ 1 ■ • 1 

m„ times 



1 — , and rj G [0, 1) and mj G Nq, j = 1, . . .,rc are defined by 



A = (2-ry_i) +mj + r h 



when j —I G K, 
otherwise. 



(4.11) 
(4.12) 



If j G K, the 2 x 2-block matrix Bj 



is left out. Note that rj = exactly when j G K. 



The following result characterizes prime spectral tetris tight frames. 

Theorem 4.11. Let N,M G N such that X :=M/N >2be given. Suppose <3> is a spectral tetris 
FUNTF ofM vectors in Mr associated with X = M/N > 2. Then, either 4> is prime, or 4> factors 
as 

L 

$ = 4»i U |J E N , 



i=i 



where ^> i is a prime FUNTF ofM — LN vectors and £V = { e j ) - = j is the standard orthonormal 
basis for Mr 7 . Moreover, 4> is divisible exactly when the eigenvalue X of the frame operator of 
$ satisfy 



jX-[(j-l)X\ >3 forall\<j< 
In particular, 4> is divisible if X > 3 — jj, i.e., M > 3N — 1. 



N 



gcd (N,M)' 



(4.13) 
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Proof. Let {ej}j =1 denote the standard orthonormal basis of Mr, and let 4> be a spectral tetris 
frame. We will use the notation introduced above. We claim that 4> is divisible if and only if 
ej E 4> for all j = 1, . . . ,N. If ej E 4> for all j = l,...,N, then 4>\ {ej} N j=x and {e/}y =1 are tight 
frames hence <I> is divisible. This shows one direction of the claim. Now, assume that e J() ^ $ 
for some jo = 1, . . . ,N. In other words, we assume that m, = 0. We consider the two 2x2 
blocks Bj -i and Bj . Let io + l E { 1 , . . . , M] for / = 0, 1 , 2, 3 be the indices of the columns 
of $ associated with B/ _i and B/ . Let <p 70 denote the j'oth row of <I>. Assume towards a 
contradiction that 4> = 4>i U<3>2 is divisible, where 4>i = {<Pi} ieIl and 4>i = {<p ; } ;e/2 are tight 
frames. The common support of the rows (p J0 and (pi° +l is {/o + 2,z'o + 3}. Therefore, owing 
to orthogonality requirements of the rows, the frame vectors <p ;o +2 and <p,- +3 need to belong to 
the same divisor, say 4>i. This in turn means that <p ;() , <p, +i E 4>2 since 4>2 otherwise cannot 
span M N . The square norm of the j'oth row of 4>i is fly o + <n o = fj , and the square norm of 
the (jo+ l)th row of 4>i is at least ^ () + (—bj ) 2 = 2 — ry . Since $i has equal row norm, this 
implies that rj > 2 — rj , that is, rj > 1 which contradicts rj E [0, 1). Hence, 4> is prime. This 
completes the proof of the claim. 

In [|23l it is shown that ej E 4> for all j = 1, . . . , N if and only if 

(j -k n )X-[(j -l-k n )X\>3 (4.14) 

for any 70 such that k n + 1 < jo < k n+ \ for all n = 0, 1,. . . ,gcd{./V,.M}. However, if (14.141) 
is satisfied for one n = 0, 1,. . . ,gcd{A^,M}, then it is satisfied for all n = 0, 1, .. . ,gcd {iV,M}. 
Hence, e ; G for all j = 1 , . . . , N if and only if 

JoA-L(j'o-l)AJ >3 

for all 70 satisfying l=k + l< j < k\ = N/ gcd{N,M}. □ 

The prime FUNTF 4>i from Theorem 14.1 H is actually the output of the spectral tetris algo- 
rithm with M — LN vectors. The frame bound of 4>i is X\ = = A — L. Note the X\ might 
be strictly smaller than 2, e.g., for M = 1 1 and N = 4 considered in the introduction, one finds 
that L = 1 thus X\ =7/4. This should be compared to the fact that the spectral tetris algorithm 
is only guaranteed to work when M > IN. It is not difficult to see that the construction also 
works when M > IN — 1 , but it is in general difficult to see for which M and the construction 
still works when N < M < IN — 1 . The next result characterizes when this is indeed possible. 
We denote the output of the spectral tetris algorithm STF(/Y,M) for all N,M E N. If the spectral 
tetris does not work for the given /V and M, we set STF(7Y,M) = 0. 

Corollary 4.12. Let N,M EN be given such that N < M < IN. Define M = M + N and X = 
M/N > 2. Then STF(/Y,M) is a FUNTF if and only if X,N and M satisfies fi4J3\> . 
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